Shell-model description of monopole shift in neutron-rich Cu 
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Variations in the nuclear mean-field, in neutron-rich nuclei, are investigated within the framework 
of the nuclear shell model. The change is identified to originate mainly from the monopole part of 
the effective two-body proton-neutron interaction. Applications for the low-lying states in odd-yl 
Cu nuclei are presented. We compare the results using both schematic and realistic forces. We also 
compare the monopole shifts with the results obtained from large-scale shell-model calculations, 
using the same realistic interaction, in order to study two-body correlations beyond the proton 
mean-field variations. 

PACS numbers: 21.10.-k, 21.30.-x,21.60.Cs 

o : 

, The nuclear shell model is based on the assumption that each of the nuclcons in an atomic nucleus moves indepen- 
PsJ ■ dently in an average field created by the other nucleons. One of the major successes of the model, even in its simplest 
version when the average potential is approximated by the the harmonic oscillator potential with the centrifugal term 
^ . (I ■ I) and spin-orbit interaction {I ■ s), was the reproduction of the so-called magic numbers, i.e. the numbers of pro- 
pLn " tons and neutrons which give additional stability to a nuclear system (2,8,20,28,50,82,126). Inclusion of correlations 
, beyond the mean-field approximation through the introduction of the residual two-bo^' interaction provides us with 

an excellent description of atomic nuclei in the vicinity of the valley of /3-stability 0, Q| . 
[ The shell structure of nuclei with large proton or neutron excess differs considerably from that of the stable 
T-H ' nuclei. The important role in the mechanism of changes is played by the proton-neutron interaction as has long 
J> I been recognized 0j 3 ■ It was demonstrated by the mean-field calculations @ that the neutrons in very neutron-rich 
nuclei experience an average field with a much more diffused boundary due to the close-lying continuum and increased 
importance of the pairing interaction. As a result, the shell structure resembles that of a harmonic oscillator with the 
spin-orbit term only. 

^— ^ ' Within the shell-model concept, the description of the nuclear mean field is conventionally obtained by considering 
, the so-called monopole Hamiltonian constructed from the centroids of the two-body interaction 6, 7, 8, 9]. The 
• eigenvalues of this Hamiltonian were shown Q to provide the average energies of the specific spherical configurations. 
' They are usually referred to as effective single-particle energies, and should reproduce the energies of single-proton 
I (single-neutron) states in odd-A nuclei with Z (TV) equal to a magic number plus or minus one proton (neutron). 
^Jh ' Recent shell-model calculations of neutron-rich nuclei with masses below A « 40 studied the mechanisms of the 
^ : changes in the shell-structure of the neutron mean-field [lol[T]| . It was shown that an increasing N/Z ratio, certain 
neutron orbitals move up or down thus deteriorating the usual shell closures and creating new ones. For example, the 
well-pronounced shell gap at iV = 20 in nuclei near ^"Ca, disappears for nuclei with Z < 12, while V = 16 instead 
appears as a new magic number. This is a consequence of the crucial role played by the attractive proton-neutron 
effective nucleon-nucleon interaction (especially, in T = channel) in which spin-orbit partners are involved |3 • 
It can be expected that the proton-neutron interaction will influence in a major way the effective single-proton 
' energies of the proton mean-field. It is well-known that the ordering and spacing of the proton mean-field orbitals in 
neutron-deficient or neutron-rich isotopes is not the same compared to /3-stable nuclei, and may exhibit pronounced 
differences. 

In this article, we intend to investigate the evolution of the proton mean-field from nuclei, situated around the valley 
of /3-stability, towards very neutron-rich systems. In particular, we consider nuclei near semi-magic ones, namely, those 
which differ from a semi-magic system by an extra proton particle (hole). Typical examples are the F, Sc, Cu, Sb and 
Bi isotopes. 

The advantage of this choice is that the low-lying states in these nuclei usually carry a major part of the proton 
single-particle strength and thus a direct comparison between theory and experiment becomes possible. These nuclei 
can be rather well described by an inert core with one extra proton particle (hole) and a number of valence neutrons 
outside the core. We then estimate the effective single-particle energies (e^) in the proton mean- field as a bare single- 
particle energy (e-,) plus the average proton-neutron part of the interaction 12, 13, 14] . In this approximation, the 
evolution of the single-particle states as a function of neutron number will depend on the specific proton-neutron 
effective interactions ((iu-ii/l^ljirji/)) and the pairing properties amongst the neutrons. 

The existence of the single-proton (single-neutron) shifts is well known experimentall y in a series of isotopes (iso- 
tones) ^3 . Recently a lot of new data has been obtained on odd- A F [161 ll^ IT^ , Cu and Bi [20l l2ll | isotopes, 
up to very neutron-rich species. 

In the present study we apply this approach to neutron-rich Cu-isotopes. We choose two types of the effective 
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interaction, a schematic zero-range force and a reahstic one, based on the G-matrix |22l |23|. in order to test the 
sensitivity to the use of various (schematic versus reahstic) forces. Finally, a comparison of the monopole shifts with 
the results of the large-scale shell-model diagonalization is carried out. 



The interest in neutron-rich nuclei around ^^Ni and in particular the questioned rnagicity of Ni has recently 
motivated the global study of low-energy nuclear structure for both odd-A Cu nuclei 13| at the LISOL facility in 
Louvain-la-Neuve and odd-odd mass Cu at ISOLDE, Cern [23| up to very neutron-rich species. Figure 1 contains the 
selected states of odd-^ Cu with spin and parity quantum numbers of the valence proton oscillator shell above Z — 28: 
2p3/2i 1/5/27 2pi/2i 159/2- The ground state of odd- A Cu- nuclei up to '''^Cu is 3/2 ~ conform with the occupation of 
the 2p3/2 oscillator orbital. As seen from Fig. 1, the 5/2^ state starts to drop drastically in energy for nuclei with 
A > 69 and might be expected to become the ground state in heavy odd-^ Cu. 

Since we are interested in the variation of the proton single-particle energies, it is consistent to compare our 
theoretical results with the experimental centroids of the 2^3/2, I/5/21 and I59/2 single-proton states, defined as 

E{nlj) = (1) 

where are spectroscopic factors for a proton transfer (from (nlj) single-particle state), and E^, - is the experimental 
excitation energy of the fcth state with J = j (both taken from [H, l^l 127, ^S. and Ref. jl^). 

Unfortunately, the data on spectroscopic factors is rather scare and differs significantly from one experiment to 
another. In Table I we summarize the available spectroscopic strength for the states of interest in ^^~^^Cu (see 
captions of the tables for references). The experimental centroids are plotted in Fig. 2(a). For those states (isotopes) 
for which no experimental data exists, only the lowest states of a given spin is indicated. As seen from the table, the 
sum rule (for the transfer to both, Tf = Ti ± ^ states) is far from being accomplished, especially, for 9/2+ states. 
Thus, one should be cautious with the trends of experimental "centroids" for this particular state. 



The low-lying states of a nucleus with one proton (neutron) outside the doubly magic core can be considered as 
mainly single-particle states and they usually serve to obtain information on the single-particle energies used in 
the shell-model. Here j stands for the quantum numbers (n, l,j, m) of a given orbital. 

The interaction between the valence proton and the filled neutron orbital results in a change of the single-particle 
energy that can be expressed as 0, |^ 



fix = + n^Ej^j. > (2) 
where Ui, is the number of valence neutrons and Ej^j^ is the average matrix element of the proton-neutron interaction: 



with V standing for the effective two-body interaction, and J being the total angular momentum of a two-body 
configuration. 

Taking into account that the model space for neutrons in general consists of several single-particle orbitals and that 
the pairing force spreads neutrons over these valence orbitals, one arrives at the following expression (12i.|[T3LIT^ : 

= + E ^^^^^ + l)vl , (4) 

where v"^^ stands for the occupation probability of the orbital ji, and the summation is performed over all valence 
neutron orbitals. 

From cq. ^ it is obvious that the proton-neutron interaction causes the proton single-particle states to change, 
depending on its character (attraction or repulsion) and magnitude. The larger the overlap between the proton and 
neutron radial wave functions and the higher ji, of a given neutron orbital is, the stronger the variations are. 

Remark that only the monopole component from the multipole representation of the effective force contributes to 
the average matrix element ^ [31 • This is why the effect is referred to as the monopole shift. 

The basic ingredients required in I^J are: (i) the shell-model space; (ii) the effective proton-neutron interaction, 
and (iii) neutron occupation probabilities. 
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For the description of the Cu- nuclei we have chosen the shell-model space consisting of four valence orbitals, 2^3/2, 
I/5/2J '^Pi/2j E^nd 1(79/2, both for neutrons and protons, outside ^^Ni-core. This model space contains all necessary 
degrees of freedom to describe low energy nuclear structure. 

In the present work we have restricted ourselves to two types of the effective interaction. First, we consider the 
zero-range ^-interaction with the spin-exchange term [3j 

= -Vo(5(7V - rv)(l - a + aa^CTj,) . (5) 

The two parameters, the overall strength Vq and the strength of the exchange term a are to be adjusted to the 
experimental data in odd-odd nuclei from the mass region one is studying. 

The geometrical properties of the zero-range interaction, even including spin-exchange contributions ((SJ, allow to 
obtain a very simple analytical expression for the monopole shift: 

where F^^ ^ ^ are Slater integrals. In order to derive this result, we had to make use of a peculiar property of the 
Wigner 3j-symbols, as described in Appendix A. 

We stress the fact that since the average proton-neutron matrix element from the (5-interaction entering jnj does 
not depend on j^^, the relative shift for spin-orbit partners (see e.g. the result for 2pi/2 and 2^3/2 orbitals in Fig. 2(b)) 
stays constant and remains equal to the difference of the original single-particle energies {ej^ — tji ). 

For the Cu- isotopes, we have determined the parameters of © to be Vb = 400 MeV.fm'^ and a = 0.1. These values 
suit the systematics known from the description of heavier isotopes '14^ and they result in two-body matrix elements 
of the same order as those of the realistic interaction which will be described below. The chosen interaction also gives 
a reasonable, although slightly stretched, spectrum of ^^Cu [33l |. 

The second type of interaction used here is a realistic one, obtained from the G-matrix through a Brueckner-Hartree- 
Fock procedure j22i] and modified further for a monopole correction 23] following the recipe of Ref. 7] to account 
for the saturation properties. The interaction works very well in the description of Ni and Cu-isotopes, reproducing 
many known spectroscopic properties |23| |. 

The neutron occupation probabilities Vj^ have been obtained from the solution of the BCS equations using a pairing 
Hamiltonian |3^: 

H — —G ^ a],„a] _^aj'm'aj',_m' . (7) 

jm.j' va' 

The strength of the pairing interaction was chosen to reproduce the pairing gaps in neighboring even-even Ni- 
isotopes |3J| resulting in a dependence G = 2?,/A MeV. 



In Fig. 2 we compare the experimental energy centroids (Fig. 2(a)) with the results of the calculations according 
to eq. (0J with two types of the interaction (Fig. 2(b,c)). 

Figure 2(b) contains the results obtained using the (5- force. For ^''Cu, the lowest states are given by the empirical 
values tj^. These values are modified by the two-body interaction. Effective single-proton energies for the I/5/2 and 
1(79/2 states first move up, when neutrons are filling the pf orbitals, and then (from ^^Cu) go down relative to the 
2p3/2 state as the neutron I59/2 orbital starts to be filled. This is governed by the relative magnitude of the centroids 
of the interaction. In particular, the drastic lowering of the I/5/2 and 159/2 orbitals stems from more attractive 
-^1/5/2139/2 ^n'i -^139/2 is9/2 values relative to -E'2p3/2 199/2- The energy of the 2pi/2 orbital stays unchanged with respect 
to the 2p3/2 ground state as has been mentioned before. At A = 79 we observe a crossing between the effective 2pi/2 
and 1/5/2 states, however, the chosen interaction does not reproduce the expected crossing of the effective 2^3/2 and 
1/5/2 orbitals. 

In Fig. 2(c), we present the results obtained with a realistic effective interaction described above. Only diagonal 
proton- neutron matrix elements contribute to the expression Q . Though the overall trend obtained from both types 
of the interaction remains the same, as follows from the similar character of the matrix elements, there are certain 
distinctions. First, the energy of the 2pi/2 orbital varies relative to the 2^3/2 ground state orbital (see Fig. 2(c)). 
Secondly, the proton I/5/2 orbital decreases steeper as neutrons are added in the lgg/2 orbital and for ^^Cu the 
inversion of the I/5/2 and 2^3/2 orbitals is predicted, compared to the use of the (5-interaction. 

We point out that the lack of agreement between the variation of the effective 1(79/2 single-proton state and the 
experimental 9/2+ level is due to the fact that the single-particle strength of the 139/2 state lies higher up and the 
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lowest state which is plotted in Fig. 2(a), due to the lack of experimental data, contains only a small fraction of that 
(see Table I). 

In order to get an insight into additional correlations that cause a variation of the single-particle centroids of 2^3/2, 
1/5/21 2pi/2 and lgg/2 proton states beyond the monopole shift, we have carried out large-scale shell-model calculations 
using the same realistic effective interaction as described above. The model space consists of the (I/5/22P3/22P1/2I59/2) 
orbitals outside the ^^Ni-core, truncated in case of ^^Cu, ^^'^^Cu and ^^Cu (up to 3, 4 and 8 particles, respectively, 
were allowed to occupy the I59/2 orbital). The calculations were done using the shell-model code Oxbash f35|. The 
single-particle centroids for the 2p^/2, 1/5/21 '2pi/2 and I59/2 states obtained by taking into account up to about 400- 
800 excited states weighted by their spectroscopic factors (in order to exhaust the sum rule for Tf = Ti — ^ transfer 
up to 98% or to get a reasonable saturation of the total value for truncated spaces) arc shown in Fig. 2(d). The wave 
functions of the corresponding Ni-isotopes have been calculated without any restrictions of the model space. 

A remarkable result is that the centroids shown in Fig. 2(d) are rather close to those obtained by the analytical 
approach Q) as shown in Fig. 2(c) (the wiggles in the behavior of the shell-model l(7g/2-centroid in Fig. 2(d) for 
65-69(^^ are related most probably to the truncations of the model space which is otherwise too large to be treated 
numerically). The monopole energy shifts incorporates the effects of a certain class of residual interactions (self-energy 
correction of the proton single-particle energy because of proton-neutron interactions and neutron pairing) on the bare 
proton single-particle energy within the model space, as expressed by eq. Q. 

A comparison of these results thus illuminates the question what type of correlations is important in the calculation 
of the energy shifts beyond the monopole shift. A rather good agreement between the centroids of Figs. 2(c) and 2(d) 
confirms that the basic ingredients of the single-particle centroid variations are already taken into account by eq. 10} . 
This is a most interesting and valuable result. 

To illustrate the importance of the pairing effects, we plot in Fig. 2(e) the effective single-proton energies obtained 
from the monopole Hamiltonian with normal filling of the single-particle orbitals by neutrons (the approach chosen 
in Refs. S S llfl Ii3 ) ■ H seen that through the occupation probabilities, the pairing force (Fig. 2(c)) smoothes 
significantly the single-particle behavior. Fig. 2(f) shows the same as Fig. 2(e), but with a different order of filling of 
neutrons orbitals (neutron 2px/2 is filled before neutron 1/5/2)- The behavior of single-particle orbitals does not have 
a zigzag structure in this case, but the maximum of the centroids is shifted to ^^Cu, compared to the calculations with 
pairing effects (Fig. 2(b,c)). Let us note that physically the order of filling chosen in Fig. 2(e) is more realistic than 
that implied in Fig. 2(f), as follows from the experimental spectra of ^''''^^Ni (neutron single-particle orbital 2pi/2 is 
higher than I/5/2 in these nuclei). 

Since the lowest-lying J'^ — 1/2^, 3/2^, 5/2^ and 9/2+ states represent only a fraction of the total proton single- 
particle strength, we compare the energies of the lowest l/2~, 3/2~, 5/2^ and 9/2+ experimental states (Fig. 3(a)) 
and theoretical states obtained from the large-scale diagonalization (Fig. 3(b)). In general, we remark a very good 
agreement with the experimental data due to the quality of the effective interaction. 

By comparing the results of Figs. 2(c,d) and 3(b), it is clear that the correlations beyond the monopole approximation 
are necessary to understand a detailed fine structure of the energy variations for the lowest l/2~, 5/2~ and 9/2+ 
levels. This is particularly so for the 1/2^ and 5/2^ levels where specific kinks, related to a possible shell closure at 
N = 40, are clearly seen in the diagonalization. The general trend of the 9/2+ level, as obtained from diagonalization, 
is quite different compared to the monopole shift. These effects are obviously a result of interactions of the pure 
single-particle configurations with many of the excited configurations with the same J'" value. 

Finally, in Fig. 4 we summarize the spectroscopic factors of the 2^3/2, I/5/2, 2pi/2 and 139/2 proton stripping 
from the '^"-'^Ni-core into the lowest-lying l/2~, 3/2~, 5/2~ and 9/2+ states of "^Cu, respectively, as provided by the 
shell-model diagonalization. The spectroscopic factors have been normalized to the sum rule presented in Appendix 
B. The absolute spectroscopic factors for the lowest states of a given spin are shown in Table I in comparison with 
the available experimental data on ^^~^^Cu. The calculated spectroscopic factors are in rather good agreement for 
the lowest 3/2~ and l/2~ states (within 20%), however, the disagreement increases for 5/2~ and 9/2+ states in these 
isotopes. It is not clear, at present, what the origin of these differences is. It may be due to a lack of experimental 
data, but it may also point out certain deficiencies in the realistic interaction presently used. As for the whole series of 
Cu-isotopes, the calculations show that only the 3/2~ ground state remains more or less a pure single-particle state. 
The lowest l/2~ state is predicted to carry about half of the 2pi/2 single-particle strength. The lowest 5/2~ state 
shows much more mixing for ^^~^^Cu. This is in line with experimental observations of a few 5/2^ states in these 
isotopes even at low energies. Finally, the 9/2+ state is predicted to be a rather pure single-particle 139/2 state only 
for fight isotopes. Starting from ^''Cu, the single-proton strength almost vanishes. In order to test these predictions, 
more precise data on spectroscopic factors and level spin and parity assignments would be very helpful. 

In conclusion, it becomes clear that using both simple effective and realistic forces in order to derive the corre- 
sponding monopole energy shifts, the experimental proton single-particle centroids are, at present, only moderately 
described. This comparison is hampered by the lack of reliable and extensive experimental results on one-nucleon 
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transfer reactions. 

An interesting point is a reproduction of a lowering in the I/5/2 energy centroid, once the = 40 neutron number 
is passed, in particular using realistic forces. We also remark that correlations that define the monopole energy shift, 
as defined in eq. Q , are very much consistent with the full theoretical single-particle energy centroid resulting from a 
diagonalization of the full energy matrix, encompassing all possible residual interactions within the large model space. 
This result is an interesting one and has to be studied in other mass regions, too. Although the energies obtained 
from the diagonalization in the large model space concerning the lowest state of each spin, parity, following rather 
well the corresponding experimental energies and exhibit interesting fine structure, the theoretical centroid, using a 
significant number of the eigenstates, becomes rather smooth again. At the same time, more experimental data on 
structure of nuclei from the region would serve a good testing ground for the present results and the interactions used. 
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Appendix A: Properties of the (5-interaction 

The diagonal matrix element of the interaction ^ between the two-body states in jj-coupling reads 

(»; J\V\m.; J) = ~Vo{2j. + 1)(2> + 1)((1 -2a- 2a(-l)'0 ^J^ J^^^ ) ^ ( ^2 1/2 -^1 ) ^ ' 
In order to evaluate the average matrix element. 



using the above two-body proton- neutron matrix elements, we have to exploit the well-known summation properties 
for 3j-symbols, 



E(2^+l)( //2 -t/2 0) =1' (10) 



L/2 -1/2 
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E(2^+1)(V2 V2 -^l)' = l- (11) 
Besides, in the calculations the particular expression, 

E(-1)'^(2^+1)(V2 -1/2 0)'' ^''^ 

appears for which we have found no standard result in the literature on angular momentum algebra (e.g., Ref. |2fl]). 
We had to be able to establish that the above expression (fT^ vanishes, as in the most general case, 

due to the peculiar summation properties of the considered 3j-symbols which we would like to mention here: 

E (-1)^(2^+1) -m 0)' = '/'' ^''^ 
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J— odd 



E(-l)^(2^+l)(i 1 ) (15) 
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Such summations on the odd or even values of the angular momentum are often encountered in quantum mechanics. 

As a result, we obtain for the average matrix element of the proton- neutron interaction the following analytical 
expression: 

E,^j. ^ ~VoF^,^i^.^„^i^il - a) , (16) 

where i i stands for the Slater integral. 

Remark, that the average matrix element does not depend on the j quantum number of the orbitals involved, as 
well as it vanishes for the spin-exchange term. 

In a similar way, we can evaluate the average particle-hole matrix element. The particle- hole matrix element relates 
to the particle-particle one via the Pandya transformation 32]: 

{mZ'\J\v\m-': J) = - Y^{2J' + 1) I J I Ow.; J\v\m.; J) ■ (17) 

Exploiting the following summation properties for 3j-symbols, 
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and making use of the peculiar relation (|12|l , one obtains for the average particle-hole matrix element the result 

Ej^j-^^yoFli^.^^^,^{l~a) . (20) 

This expression differs from that for a particle-particle average matrix element only by its sign. 

Appendix B: Sum rule for the spectroscopic factors 

The total spectroscopic strength in proton stripping reactions to the states with Tj ~ Ti — ^ is given by the following 
sum rule Is^ : 

E §TT^^^"^^'^ ^ 2TUm + 1) ^ ^ - - + ' (21) 

JijJf are the spin and Ti,Tf are the isospin quantum numbers of the initial and final states (in ^""'^Ni and '^Cu, 
respectively), C — {Ti,Tiz, \, —\\Tf = Ti — \,Tfz = Tiz — \) is the Clebsch-Gordan coefficient, (n(j)} is the average 
number of proton or neutron holes in the j-orbital in the initial state (ground state of "*~^Ni in our case). 
Since in the case considered here, Ji = 0, Jf = j, the sum rule (|21|l reduces to 

E Sfinimjf + 1) = (,V(j))t. - ^f^M))T. , (22) 

as indicated in Table I where the average number of neutron holes in the ground state of Ni-isotopes are taken from 
the shell-model diagonalization. 
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TABLE I: Available experimental data on spectroscopic factors of the lowest 3/2", 5/2~, l/2~ and 9/2+ states in 59.61>63,65q^ 
(taken from taken from Ni(d,n) or Ni(^He,d) proton-stripping reaction studies |23. l2ft I27ll28. 29, 30]). The sum rule refers to 
the proton transfer to both, = Ti ± i states, while the shell-model sum rule is given for the states with Tf — Ti — ^ only. 
The calculated spectroscopic factors for the lowest states of a given spin are indicated as well. 
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FIG. 1: Experimental low-lying states in odd-yl Cu isotopes [l9ll25|. 
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(d) centroids from_diagonalization 
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FIG. 2: (a) Experimental centroids when data is available (see Table I) or lowest states of a given spin in odd-A Cu isotopes; 
(b,c) theoretical evolution of the single-particle states in odd-A Cu isotopes duo to the nionopole part of the proton-neutron 
(5-interaction and realistic interaction, respectively; (d) shell-model centroids derived from the exact diagonalization, using a 
realistic interaction; (e,f) the same as in (c), but implying normal filling of the neutron orbitals (two possible ways of filling). 
The arrows in figure (b) below the plot indicate schematically the normal filling of neutron single-particle orbitals in Cu-isotopes. 
See text for details. 



11 



> 

(U 



LU 



(a) experiment 



9/2 



9/2 



1/2 



3/2 



f-t-y-f-t-t-f-t-^3/2- 



55 57 59 61 63 65 67 69 71 73 75 77 

A 



4- 



> 
0) 



LU 



(b) exact diagonalization 

9/2' 



K 



5/2 



1/2" 

p 



1/2>o «?;o-o 
3/2"b-0'^° ^A_«^'^ 

f-y-f-f-f-?-f-f-t-y-r*f5 /2 



57 59 61 63 65 67 



71 73 75 77 79 



FIG. 3: The lowest 3/2 , 5/2 , 1/2 and 9/2+ states in odd-A Cu isotopes: (a) experimental values; (b) as obtained from the 
shell-model diagonalization with the realistic interaction. See text for details. 
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FIG. 4: Proton stripping spectroscopic factors into the lowest 3/2^, 5/2^, l/2~ and 9/2+ states, respectively, in odd-A 
Cu isotopes obtained from the shell-model diagonalization with the realistic interaction. The values are normalized to the 
shell-model sum rule. See text for details. 



